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ing in the frame of the simple quark-parton model. 
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Introduction. 



The increasing amount of evidence for i/^^^i'^ oscillations entails that some 
experiments being prepared should actually 'see' the r produced by C.C. 
Vt interaction. Accordingly, papers have appeared of late with the purpose 
of giving the experimental groups the necessary tools to prepare their Monte- 
Carlo and the analysis of their future data. One of the themes of these papers 
is of course the calculation of the r production cross-section, both differential 
and integrated, and that of the r polarization which is necessary to predict 
the angular distribution of the decay products into various channels, see e.g. 
[1][2][3] 

However, not all these papers give simple recipes which can easily be imple- 
mented in simulation programs. We therefore thought it useful to write this 
note which describes the way we calculated the cross-section and polarisation 
and implemented them in the simulation which served many years ago to pre- 
pare the proposal for the NOMAD experiment [8]. The following is essentially 
an improvement over an unpublished internal note [9] 

We shall only treat the deep inelastic scattering (DIS) case in the quark- 
parton model. This makes things particularly simple since quarks are then 
described by Dirac spinors. However, the procedure can be straightforwardly 
extended to quasi-clastic or resonance production, through the introduction of 
form factors. These complicate the expression for the hadronic current, but the 
main result, viz the extraction of the r polarization 4- vector can be done by 
mimicking what follows. 

Formulae for C.C. neutrino interactions without neglect of the lepton mass 
were first published in [5] (see also [6]). However, they were written in terms of 
the scaling functions i^/s and no connection was given with the basic quark/anti- 
quark densities given by the available structure functions library packages. 
Moreover, [5] gives a 'polarized' part of the cross-section which is not very 
useful for Monte-Carlo, whereas the polarization vector to be used in decay 
matrix elements is not explicitly written. It could be extracted from the given 
formulae but the result would only be an average with respect to the vari- 
ous elementary nuclear constituents. In using Monte-Carlo of the Lund type, 
where proper simulation of the hadronic final state requires knowledge of the 
type of qTiark/anti-quark which has been hit, it is more appropriate to derive a 
polarization for each case. 

In the spirit of the impulse approximation, wc shall therefore start from the 
elementary interactions and derive what should represent the observable cross- 
section and polarization by summing and/or averaging. We also make explicit 
the derivation of the kinematical limits in x and y which are only stated in the 
above mentionned papers and can be used as training for the students. 

Notations and basic formulae are listed in section 1; in section 2 we derive the 
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kinematical bounds for neutrino-nucleon C.C. scattering with a massive lepton 
in the final state. Going over to the quark-parton model, the elementary sub- 
processes cross-sections and polarization 4-vectors for the lepton are calculated 
in section 3. Section 4 sums up these results with a neutrino-nucleon cross- 
section formula and an average lepton polarization. 



1 Notations 



The reaction to be described is i^r ^ ^ t X. 
Four momenta : k = PviP = Pn,Ptj Q = Pv — Pt 
Masses : rur for r , M for target nucleon, W for X. 
Kinematical variables: 

C.M.S. energy squared: s = AI'^ + 2ME with E the neutrino energy in the 
target rest frame. 

Four- momentum transfer squared: Q'^ = —q^ 
Squared mass of final hadronic system: W"^ = (g 
Bjorken's variables x: x = /2p.q y =p.q/p.k 
C.M.S. scattering angle: 6 
Useful relations: 



Q2 = xy{s - M^) = -ml + 2Ei-"'- {E^-""- - f^-^ cosQ) (1) 

W'^ = s + Tnl-2E''^'^-^s (2) 

_W^-M^ + Q'^ _ W^-M^ 

^~ S-M2 ~ (l-ar)(s-M2) ^' 



2 u — N ^ I — X kinematics with non zero lepton 
mass. 

From (3) and > one finds a; < 1 as usual. Replacing W'^ taken from 
(2) into (4) leads to the following expression for £;^.m. . 

pern. _ ml + js - M^){1 - y + xy) 

- 

On the other hand, the second of relations (1) and \cos6\ < 1 yield: 

|xy(s - M2) +ml- 2E^-'^-E^-"'-\ < 2p^-™-^^-™- 

Now define Sj- = jzm^^ = s-m^ ' E^-"^- = to rewrite this: 

rpc.m. cm. 

\xy + 5r-^\<^ 
\/s -s/s 
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or letting h xy + Sr and squaring: 

Tpc.m. ™2 

h'-2h^ + ^<0 
\/s s 

Using (5) for E'^ "*- and the definitions of (5t-,(5jv, this can be transformed to: 

h^5N -h{l-y)+5r <Q (6) 

or, by re-expressing y as function of h: 

(1 + x6N)h^ -{x + Sr)h + xSr<0 (7) 

Both inequahties (6) and (7) lead to the kinematical limits in terms of x and y, 
albeit in slightly different forms. Solving the first one for h will give limits on 
X clS cl function of y and conversely for the second. This later inequality leads 
to the limits quoted in [5] as follows: (7) is possible only if: A = (x — Sr)^ — 
ix'^drSN > 0, i.e.: 

\x-6r\> 2x^/6r5N (8) 

And the limits on y for given x are then found from the roots of the trinomial 
in (7): 

2{1 + xSn)x 2{l + xdN)x 

Clearly, one must have x > 5r so that (8) is to be understood as : 

1-2VSJ^ 1-2^ 1-^ ^ ^ 

Note that the upper limit for x derived from (6) is irrelevant (above 1) but that 
as soon as one demands a minimum value for W above M (DIS should mean 
W > Ma at least), then relation (4) entails an upper limit on x as a function 
of y 



3 Elementary interactions. 

To derive formulae for the differential cross-section and polarization, we use the 
simplest quark-parton model and then sum over partonic contributions. 

3.1 Matrix elements. 

The kinematics is pi, +pq = Pr +Pq' ^ and use shall be made of the r polarization 

four-vector S 

The reaction amplitude is: 

^iZ(r)7"(l - j')u{u)u{q')j^{l - j')u{q)B.W. 

^we will use this notation even for anti-quarks and anti-neutrinos except in the 
final x-section formulae 
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where 

BW = ^ 

Ml^-q^ -iVwMw 

is a W's propagator correction to the pure Fermi amphtude, where terms arising 
from the longitudinal (^r^) part have been neglected. 
Squaring and introducing the density matrices one gets: 

^Tr(p,7"(l - 7')P.7^(1 - 7'))rr(p,,7a(l - l^)p<ilp{-^ - l')) 

Here pr = 5(7 -p,- + "ir)(l + 7^7 ■ 5') but the other density matrices are summed 
(averaged) over polarizations and is understood here and in what fol- 

lows. 

The first trace is the so-called leptonic tensor, which, with pi, = j ■ k {p^, = k 
here, for ease of notation) reads explicitly: 

where L k — rrirS 

For anti-neutrinos, L should read k + rrirS and the sign of the antisymmetric 
part should be reversed. 

Taking the trace of the quark tensor and contracting, one finds: 

|rp = 32Gpp^ ■ pq{pr - rrirS) ■ Pqi (11) 

We shall use the S dependence of this result later. For the moment, we sum it 
over polarizations to get the transition probability: 



v ' PqPr ' Pq' 

For an anti-quark target, one gets instead of (11): 

32G|p^ • p^{pr - rrirS) ■ pg (12) 

so that the roles of the initial and final partons are permuted. Finally for an 
anti-neutrino going over into r+ one finds: 

32G|npj7 • Pg' {pr + rrirS) ■ pg (13) 

and 

32GIpjj ■ pgipr + rrirS) ■ p^ (14) 
for a quark and anti-quark target respectively. 
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3.2 Cross-sections. 



The differential cross-section for each of the elementary processes considered 
above is: 

d'^=p\T\ (27r) 6 {Pq' +Pr-Pq-Py)^^^ 

with F the Moller flux factor and \T\ is the matrix element computed above. 
Standard manipulations and integration with respect to the r azimuthal angle 
reduce this expression to 

da = — — iTPdy 

with y defined in section 1 as the leptonic fractional energy loss in the target 
nucleon rest-frame. The flux factor is simply F = Ap^ ■ pq. 



3.3 Polarization 



Here shall be found the main simplification of our presentation w.r.t. others. 
Setting up a polarization basis is useful for testing symmetries or conservation 
laws in the production process, but is of no use for simulations which only 
require the r polarization four-vector in any specified frame, to be fed, for ex- 
ample, in a decay matrix element. 

The calculation is most simply done by following the reasoning of ([4]): the 
squared amplitudes written above are proportional to the probabilities for find- 
ing the final in a given state of polarization described by the four-vector S. 
If the true polarization from the production process is 5/, then those proba- 
bilities are equally found by projecting the true density matrix 1 -|- 7^7 • Sf on 
1 -l- 7^7 • S which represents the given state for which we want the probability. 
In other words, \T\^ oc Tr{l + 7^7 • >S')(1 -|- 7^7 • Sf) oc 1 — 5 • 5/ 
Hence we find for the v — q case: 

^ Pt ■ Pq' 

Since Sf Pr = Owe cannot simply invoke the arbitrariness of S to cross it away 
but we must allow for a term proportional to Pr'- 

nirPq' 

Sf = — + Xpr 

Pt ■ Pq' 



A is now determined by projecting this equality onto Pr and found to be — ^ 
so that the final result reads: 

_ nirPq' Pr 



Pr ■ Pq' rUr 



(15) 



In the T rest-frame we see that the time component 5° = as it should be and 
that the space part (the polarization in its usual sense) is: 



Eqi 
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For the three other cases one finds: 



u — q 


u-q 


V — q 









It is important to note that what has been done here with the simple quark 
tensor can be adapted to more comphcated forms of the hadronic current and 
of the tensor built from it. For example, if r is produced through a quasi- 
elastic C.C. interaction, the matrix element of the hadronic current between 
the neutron and the proton states is described by six form factors (two of which 
are zero in this case) but exactly the same procedure can be used to identify 
Sf in terms of 4-momenta, kinematical invariants and these form factors. 



4 Interactions with Nucleons 
4.1 - Nucleon cross-section 



The elementary cross-sections written in 3.2 read explicitly 

2(^2 

da{v -q) = —{Pv -PqPr- Pq')dy 

T^Pv -Pq 

da{v -q) = — {py -P^Pt- Pq)dy 

TfPu-Pq * 
2(^2 

da{v -q) = —{Pu ■ Pq' Pr+ ■ Pq)dy 

VrpF -Pq 
2(^2 

daiV -q) = — {pu ■PqPT+ ■ P^)dy 

T^Pv-Pq ^ 

Let ^ be the nucleon momentum fraction carried by the struck (anti-) quark: 
Pq = ^p. Energy-momentum conservation says that: 

Pu + ^P = PT+Pq' (16) 

or q + ^p = pg, (17) 

by squaring (17) we get: 

2{^-x)p-q + eM^ = ml, (18) 

In Bjorken's limit, q^ — > oo,p ■ q ^ oo (18) shows that ^ = x 
Squaring now (16), combining with (18) and using ^ = x yields: 



■ pqi = x{s — M^) — 
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Therefore : 

9/^2 r 

da{v-q)=da{V-q) = ^MEx{l-^)dy (19) 

Taking the scalar product of (16) hy p^, using again ^ = x and the definition 
of q yields: 

2pq, ■ = x{s - M2)(1 - ^ - y) 



On the other hand: 
and 

Hence: 



2pq ■ pr = xil - y){s - M^) 
F = Ap^ ■ Pq = 2x{s - M'^) 



2(72 r 

da{u -q)= da{V - q) = ^MEx{l - — - y)(l - y)dy (20) 

TT X 

Formulae (19) and (20) solve the cross-section question. Multiplying them by 
the appropriate quark/anti-quark distribution functions, summing over flavors 
and re-introducing the W boson propagator factor, one finds for a neutrino 
beam: 



da _ 2G\EMx 



dx dy TT 



(1 - ^)Q(x, g2) + (1 _ ^ _ _ g2 



2\\2 



\B.W.{Q')\ 



(21) 

where Q and Q are the appropriate mixtures of (Q.C.D.- evolved) quark dis- 
tribution functions for the nucleon which is hit. 

For an anti-neutrino beam, the coefficients of the quark and anti-quark distri- 
bution functions must evidently be exchanged and the distributions themselves 
adequately modified for the given nucleon, in order to take into account charge 
conservation at the constituent level. 



Expressions (19)(20) and/or (21) are directly usable with standard quark 
distribution function libraries; from that point of vuc, they arc more practical 
than expressions found in the quoted articles where nucleons are described by 
structure functions or their scaling limits. In this same limit, the connection 
between the two descriptions is easily made along the lines of what can be 
found, for example, in ([7]) for the electromagnetic interactions: one merely has 
to identify the general hadronic tensor with the sum of quark and antiquark 
tensors decomposed on the same basis; the various relations F2 = 2xFi, F4 = 0, 

= Fi follow and after simplifications, formula (21) is retrieved. 



4.2 T polarization in j/-nucleon scattering 



The two terms in (21) are to be interpreted as the relative probabilities for scat- 
tering from a quark and an anti-quark in the target nucleon. The r polarization 
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vector is therefore the weighted average of the two relevant Vf of section 3.3. 



Let us call Pq = (1 - and Pq = {I - y){l - y - 

The coefficients of Pq and Pq were found in section 3.3 to be: ^ and 
Rewriting the denominators in Lorcnz-invariant form as: 



^, , ^5 x(l - ^) and £'a — — ^ — 

we find the average r polarization in neutrino scattering: 



x{l-y) 



2mr 



and in anti-neutrino scattering: 
2mr 



{p+ ^)Q(x,g2) -"-L- y)Qix,Q^] 

X X 



M2 



Sr 



(Pq + Pq) 



-1 



with the same proviso as in the preceding subsection as to the contents of Q 
and Q 

In these expressions, all 3-vectors are expressed in the r rest-frame. Should 
one need Sf in another frame, one could either boost the four vector {0,V) 
or perform the same averaging directly with the Sf vectors extracted from the 
transition probabilities as done in section (3.3). It is clear that the relative 
probabilities are made of Lorentz invariants only and therefore valid in any 
frame, as it must be physically. 

Let us remark, however, that for a Lund-type Monte-Carlo where the struck 
quark is identified, no averaging takes place at the event level and one uses 
directly the results of (3.3). 



5 Conclusion 



We have given lowest order but simple formulae for cross-section and polar- 
ization of T 's produced by i^r 's charged current interactions valid in the DIS 
regime. We leave it to the interested reader to use her/his favorite quark distri- 
bution functions package to calculate and study cross-sections and polarizations 
with the help of these formulae through Monte-Carlo or numerical integration. 
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